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Abstract. We present a new formalism for the description of multiphoton plasmon 
excitation processes in jellium clusters. By using our method, we demonstrate that, in 
addition to dipole plasmon excitations, the multipole plasmons (quadrupole, octupolc, 
etc) can be excited in a cluster by multiphoton absorption processes, which results 
in a significant difference between plasmon resonance profiles in the cross sections for 
multiphoton as compared to single-photon absorption. We calculate the cross sections 
for multiphoton absorption and analyse the balance between the surface and volume 
plasmon contributions to multipole plasmons. 



1. Introduction 

In the present paper we demonstrate that, in addition to dipole plasmon excitations 
multipole plasmons (quadrupole, octupole, etc) contribute to the multiphoton excitation 
process, which results in a significant difference of plasmon resonance profiles between 
the cross sections for multiphoton and single-photon absorption. We have developed a 
formalism from which the cross sections for multiphoton excitation can be worked out. 
The balance between multipole surface and volume plasmon contributions is analyzed. 
Our results are obtained within a theoretical model for the multiphoton excitation of a 
jellium cluster. This model is applicable to metal clusters, to fullerenes and to any type 
of cluster in which a strong derealization of valence atomic orbitals takes place. The 
theoretical formalism we have developed is not confined in its application to photons. 
It can also be used to describe any kind of higher order plasmon excitation processes, 
for example those which arise by multiple scattering of electrons within the cluster. 
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Recently, a number of papers have discussed metallic clusters |1|, |2| and fullerenes 
H in strong laser fields. The theoretical approach usually followed is to solve a time- 
dependent local density equation numerically (LDA or local-density approximation), 
and the regime most commonly studied involves intense, short laser pulses, for which 
the turn-on and turn-off properties significantly affect excitation. 

We report here on a different problem: our initial interest lies in lower laser powers, 
for which multiphoton excitation just begins to intrude, and we are interested in 
developing a formalism to describe the interaction between collective modes and a laser 
field in the multiphoton regime. For this purpose, a semiclassical model, in which the 
collective flow of charge is driven by a periodic field, is established, and we relate it to 
the multiphoton absorption cross section of the cluster, which takes account of quantum 
mechanics. Of course, one can in principle extend this treatment to consider the turn-on 
and turn-off of laser pulses, to treat the interaction numerically for various power levels 
and initial charge distributions. What we wish to point out, however, is a novel feature, 
which arises even for an infinite wavetrain interacting with a cluster (the simplest and 
most fundamental problem): multiple plasmon excitations are driven by multiphoton 
excitations. In the present paper, we explain by what mechanism this arises. 

Our approach is based on different principles from the LDA. Instead of using the 
Kohn-Sham formalism, we propose a hydrodynamic approach. In the LDA, all quantities 
are made to depend solely on charge density, and currents are subsequently made to 
appear by solving a time-dependent equation. In our approach, we start out from the 
continuity equation and the Euler equation, from which both the current flow and the 
density are obtained in a completely consistent way. Our theory is local, and does not 
include exchange. In principle, it would be possible to extend it by building in exchange 
and correlations in a similar way to the LDA. One of the benefits of our approach is that 
all the momentum transfer terms are included in the formulation, which leads to the 
presence of both volume and surface plasmon terms. If one wishes to simplify the theory, 
it is possible towards the end of the calculation to assume zero momentum transfer, in 
which case the volume plasmon terms disappear from the problem. 

Surface plasmon excitations are well known in atomic cluster physics. The 
dipole surface plasmons are responsible for the formation of giant resonances in 
photoabsorption spectra of metal clusters (see e.g. [|], ||, [| S | ; S [Ml 0' Hll)- They 
also play an important role in inelastic collisions of charged particles with metal clusters 
fjl p"5|, |16|, [T7L [TR [I9|. The role of surface plasmon excitations in inelastic electron- 



cluster scattering was thoroughly studied in [13], [L7], [18], |19[j , and it was demonstrated 
that collective excitations make a significant contribution to the electron energy loss 
spectrum (EELS) in the region of the surface plasmon resonance. With increasing 
scattering angle, plasmon excitations of higher angular momenta become more and 
more prominent. 



3 



Plasmons are characteristic of delocalised electrons, and therefore the jellium model 
provides the most appropriate starting point for a discussion of plasmon excitation 
in the multiphoton regime. This regime is particularly appropriate for the study of 
atomic clusters, which are rather fragile objects, and readily explode under very strong 
irradiation pUf 

The inelastic scattering of fast electrons on metal clusters in the range of transferred 
energies above the ionization threshold was considered in [21]. It was demonstrated that, 
in this energy range, volume plasmons dominate the contribution to the differential 
cross section, resulting in a resonance behaviour. The volume plasmon resonances 
excited in the cluster during a collision decay via the ionization process. The resonance 
frequency and the autoionization width of the volume plasmon excitations have both 
been determined in |H[ . 



The role of the polarization interaction and plasmon excitations in the process 
of electron attachment to metal clusters has also been examined both theoretically 
22| , [23| and experimentally p4 |. It was demonstrated that plasmon excitations induce 



a resonance enhancement of the electron attachment cross section. 

Our paper is organized as follows. In section ^|, we derive quantum mechanical 
expressions for the cross sections of multipole (quadrupole, octupole etC) plasmon 
excitations taking place in the multiphoton absorption regime and estimate the cross 
sections on the basis of the plasmon resonance approximation. We present and 
discuss the multiphoton absorption profiles and demonstrate a significant change of 
the multiphoton absorption profiles in the cross sections for multiphoton as compared 
to single-photon absorption. In section |3|, we establish certain connections of the cross 
sections with the variation of electron density in the cluster due to the external field 
and present the formalism for the calculation of the electron density variation in the 
cluster due to the external field, based on the use of the hydrodynamic Euler equation 
and on the equation of continuity. We apply the general formalism to the description of 
fast electron-cluster scattering and multiphoton absorption. In section |J we calculate 
the multipole moments of the system induced by the external field on the basis of the 
formalism outlined in section |^. We analyse the plasmon resonance structure of the 
induced multipole moments and conclude that it is analogous to the one arising in the 
multiphoton absorption cross sections calculated in section |2|. In section [5], we draw 
conclusions from this work. In |Appendix A| , we show how matrix elements for collective 
transitions can be calculated on the basis of the sum rules. In Appendix B , we present 
details of calculations of the angular integrals arising in the formalism outlined here. 
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2. Plasmon resonance approximation for multiphoton absorption cross 
sections 

First, we consider the cross section for multiphoton absorption in jellium clusters and 
demonstrate that multipole plasmon excitations are essential to this process. The 
discussion in this section is based on the plasmon resonance approximation, which is 
introduced below. 



2.1. Single photon absorption 

Let us start by considering the simplest example and calculate the cross section for 
single-photon absorption in the plasmon resonance approximation. 

The single-photon absorption cross section in the dipole approximation reads as 

47r 2 e 2 

0-1 = uj^2\z on \ 2 5(uj no -huj) (1) 

C n 

Here, e is the charge of electron, c is the velocity of light, h is Planck's constant, 
^no = £ n — £o is the electron excitation energy, u is the photon frequency and ez on is 
the matrix element of the z- component of the cluster dipole moment. The summation 
over n includes all final states of the excited electron, which belong to both the discrete 
and the continuous spectra. 

In the jellium picture, which works reasonably well for metal clusters and to some 
extent for fullerenes, the main contribution to the cross section (|]) arises from a small 
group of excited states or sometimes even from a single transition of frequency close 
to the classical Mie resonance frequency - also known as the frequency of the plasmon 
resonance. For a spherical metal cluster, this frequency is given by (see e.g. |10|, [L7|] and 
section 4 of this paper) 

2 _ 4nNe 2 I 

Ul ~ ' (27+7)' [ ) 

Here V = 4nR 3 /3 is the cluster volume, where R = TqA^ 1 / 3 is the cluster radius, r 
is the Wigner-Seitz radius; N is the number of delocalized electrons in a cluster, I is the 
angular momentum of the plasmon mode, m is the electron mass. Note that, by using a 
single photon of energy 1 — 4eV, one can, in practice, excite only I = 1 dipole plasmon 
oscillations in a metallic cluster. 

For nearly-spherical fullerenes C20 or Ceo, the plasmon resonance frequency is equal 
to (see 0,0) 

_ l(l + l)N 

Ul ~ {21 + 1)723 (3) 
where iV is the total number of delocalised electrons (4 electrons per atom times the 
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number of carbon atoms in the fullerene molecule) and R is the the radius of the 
fullerene. 

The plasmon resonance approximation is based on the fact that excitations in 
the vicinity of a plasmon resonance exhaust the sum rule almost completely (see 
||, H [7], [II], [TTJ), which means that the summation in the sum rule (see e.g. P^ ) 

2_^,u no \z on \ ^ m I 4 J 

need be performed only over excited states in the vicinity of the Mie resonance. 

Now, assuming a Lorentzian distribution of width Ti for the plasmon resonance 
states and replacing the delta function, 5{oj no — Too), in ([!]) by the profile (see e.g. p5| ) 



p 

8{uj no -fuo) — > p-, (5) 

2uh{{uj 1 - uj) 2 + 

one recovers the well-known expression for the single-photon absorption cross section 
(see e.g. 1, 0) 

nNe 2 T l _ AnNe 2 u 2 Y x 

ai mc ( Wl _ w )2 + 3 ~ mc (u 2 - uj 2 ) 2 + u 2 T\ 

The width Y x is due to Landau damping. Its calculation for metal clusters is 
performed, for example, in pi] . 

The cross section (|6]) reproduces correctly the appearence of the plasmon resonances 
in single-photon absorption spectra of metal clusters and fullerenes, although some 
details of the experimentally observed profiles are naturally beyond the plasmon 
resonance approximation. The discussion of these details is not the scope of the 
present work. This can only be done accurately enough on the basis of ab initio 
many-body theories. Instead, we analyse the multiphoton absorption cross sections 
on the basis of the plasmon resonance approximation and elucidate the role of multipole 
plasmon excitations in its formation, because our interest lie in establishing the physical 
mechanisms which underpin multiphoton excitation. 

2.2. Two-photon absorption 

In the dipole approximation, the two-photon absorption cross section is equal to 

5{u no - 2Hlj) (7) 



327r 3 e 4 ft 2 _ 
° 2 = ,2 u L 



u n 



nm^mo 



\y > - / : in 



■a 



fnu — uj mo + i5 



We evaluate the cross section (|7]) in the same way as for the single-photon case. The 
main contribution to the sum over the intermediate states m arises from the virtual 
dipole plasmon excitations. Therefore, one derives 



E 



nrn^m.o 



•-.'.> in - 



m fku - u mo + \5 



2 N _ \z nl > 2 



2mhu 1 ? {u - c^) 2 + T 2 /A (8) 
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Here, we have also introduced a dipole plasmon resonance width r\ and used the 
sum rule (|J) for the evaluation of the matrix elements for the dipole plasmon excitation 
|rio| 2 ~ Nh/2muJx. The remaining matrix elements z n \ in (^) describe dipole transitions 
from the dipole plasmon resonance state to other excited states. Matrix elements for 
these transitions obey the dipole selection rule. This means that the angular momentum 
of the final state can only be equal to either I = or I = 2. According to (0) and (^), 
there is no surface plasmon excitation with / = either in metal clusters or in fullerenes. 
Thus, only transitions to the states with I = 2 are of interest. 

These arguments show that, by using two photons simultaneously, one can excite 
the quadrupole plasmon resonance in a metal cluster or in a fullerene with a frequency 
given in (||D and @ respectively. When calculating the cross section (|7|) in the vicinity of 
the quadrupole plasmon resonance excitation, i.e. at 2u ~ u 2 , it is sufficient to consider 
only transitions to the resonance final state, i.e. to put Z) n |z n i| 2 ~ \ z 2i\ 2 (here and 
below we use indices 1 and 2 to designate the dipole and quadrupole plasmon resonance 
states) and to replace the delta function, 5{uj no — 2hu) by a Lorentzian distribution of 
width T 2 (see e.g. p5|| ) 

r 2 

8{u: no - 2%u) — > p-. (9) 

2tt^((^ 2 - 2uoY + ±f) 

By substituting (H) and (|]) in (|7|), one derives 

a2 ~\ mc J 2Nhu 1 ^ LU _ iUl y + El _ 2uJ ) 2 + f 

This result demonstrates that the photoabsorption profile in the two-photon process 
differs substantially from the single-photon case. The cross section (|10J) has a resonance 
at the dipole plasmon frequency and, in addition, also contains the quadrupole plasmon 
resonance at uo = uo 2 /2. 

The cross section (|K|) depends on | ^r 2 1 1 2 - The transition matrix element z 2 \ describes 
the electron transition between the dipole and quadrupole plasmon resonance states. 
This is a single electron transition rather than a collective one. Therefore, calculation 
of Z21 on the basis of the sum rule (f|) would lead to a significant overestimate of the 
value of this matrix element. Instead, one can use Heisenberg's uncertainty principle 
for the evaluation of z 2 \ [^6[ • By estimating the radial component of the momentum of 
a single electron in a dipole and quadrupole plasmon oscillatory mode as pi ~ rmviAR 
and p 2 ~ mu} 2 AR respectively, one derives 

h h 

\z 2 i\ ~ min(Azi } Az 2 ) ~ — ~~ A — (11) 

p 2 muiAR 

Here, A is a dimensionless constant, of the order of one, Az\ ~ fi/pi and Az 2 ~ h/p 2 
are the uncertainties relating to an electron in the dipole and quadrupole plasmon modes 
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\ fn\z 2 -y\ 


1 mc . 


1 2Nhu 
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respectively, AR is the width of the domain in the vicinity of the cluster surface within 
which plasmon excitations take place. In Appendix A we prove the correctness of this 
estimate and demonstrate that matrix element Z21 is equal to 

8/6X 1 / 4 % 

Z21 ~ ~3 UJ ^AR (12) 

By substituting (|TJ) in (|TDj) , we obtain the final expression for the two-photon 
absorption cross section: 

^_/ 47riVe 2 \ 2 A 2 h ^ 1 r 2 

{ mc ) 2mu 1 NAR 2 w\ ^ _ + El _ 2^)2 + V A 

where A = | (§) 1/4 ss 2.79. 

We note that the cross section (|T3|) depends explicitly on Planck's constant h, while 
the cross section (|6]) does not. The independence of (|6]) from ft is connected with the fact 
that plasmon oscillations are a purely classical effect, while the dependence of (fl"3| ) on ft 
arises from the interaction between dipole and quadrupole plasmon modes as can be seen 
from estimate (|TT| ) and the explicit expression (O). This indicates that it is meaningful 
to treat plasmon excitations classically, while the coupling of various plasmon modes 
in the multiphoton photoabsorption process must be treated beyond purely classical 
theory. 

In figure 1 we plot the cross section profiles per unit atom for single-photon (dashed 
line) and two-photon (solid line) absorption calculated according to @ and (|13D . These 
profiles do not depend on the number of atoms in the cluster. Note that the scales are 
not the same for the two curves for reasons of definition of the cross sections in the single- 
and two-photon cases, but are both given in atomic units. The peak in the single-photon 
plot gives the location of the dipole resonance. The other peak in the two-photon plot is 
the quadrupole resonance. This figure demonstrates a significant difference between the 
nature of the profiles, arising from the presence of quadrupole plasmon excitation in the 
two-photon case. In this calculation we have input r = 4.0 and I\ = uji/4, T 2 = ^2/^1 
AR = ro- The choice of these parameters can be different for different clusters, but it 
should always lead to qualitatively similar single- and two-photon absorption profiles. 
An accurate determination of the parameters is only possible on the basis of ab initio 
calculations. 



2.3. n-photon absorption 



The formalism we have developed can also be used for the evaluation of the multiphoton 
absorption cross sections for a larger number of photons. In the dipole approximation 
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Figure 1. The profiles of single-photon (dotted line) and two-photon (solid line) 
absorption calculated according to (||) and ( |l3| ) and normalised per unit atom. The 
two-photon absorption profile is scaled by a factor 1/100. The scales are not identical 
for the two curves for reasons of definition of the cross sections in the single- and 
two-photon cases, but both are given in atomic units. 

the n-photon absorption cross section has the following form: 

(2n) n+1 n\ 2 e 2n h n -\ n 

k 

Here the amplitude M n is equal to 

,„ , ,„ l ,„ 1 ((n - \)%u) - uj mn _ lQ + iS) 

z mi o (15) 



a,, — w"53|Af fc | 5(uj ko -nhu) (14) 



M k = E E -E 



((n - 2)hu - u mn _ 20 + i5) "" (Huj - u miQ + iS) 
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The plasmon resonance structure of the multiphoton absorption cross section (p"5|) 
can be analysed in a similar way to the previous treatments for the single- and 
two-photon cases. This analysis immediately leads to the important conclusion that 
plasmon resonances with larger angular momenta (octupole, etc) can be excited in the 
multiphoton absorption regime. Thus, for example, with three photons, the octupole 
plasmon resonance at u = u 3 /3 will also be excited. This analysis, however, leaves 
undefined the matrix elements for electronic transitions between various plasmon modes. 
Estimates of these can be performed either on the basis of Heisenberg's Uncertainty 
Principle or by a calculation similar to the one for Z12 (see ( |l~2l ) and (|10|) ) in the two- 
photon case, but their accurate evaluation is not trivial. 

Note that the plasmon resonance approximation allows one to analyse only the 
plasmon resonance excitations that are characterised by relatively low angular momenta, 
because electron excitations in the cluster with large angular momenta I have single- 
particle character. This follows, for instance, from the fact that with increasing I the 
wave length of the surface plasmon mode, 2ttR/1, becomes smaller than the characteristic 
wave length of the delocalised electrons at the Fermi surface, 2Trh/\ / 2mAe, where Ae is 
the characteristic electron excitation energy in the cluster. In other words, excitations 
with angular momenta comparable with the characteristic electron angular momenta of 
the ground state exhibit single-particle rather than collective character. Therefore, when 
analyzing contributions of the plasmon resonance modes to the multiphoton absorption 
cross section, one should consider only the lowest angular momenta. For example 
according to the jellium model, the maximum angular momentum of the delocalised 
electrons in the Na^ cluster is equal to 4. Therefore, only dipole and quadrupole 
collective modes can be expected in this case. With increasing cluster size the number 
of essential plasmon modes grows as R. 

3. Hydrodynamic description of collective motion of the electron density 
in a cluster 

The multipole plasmon resonances arising in the multiphoton absorption cross sections, 
should also appear in other physical characteristics of the cluster, which can be probed 
in the multiphoton absorption regime. In the situation where plasmon resonance 
excitations are the dominant contribution to the multiphoton absorption cross section, 
it is natural to seek and analyse the plasmon resonance structure of the variation of 
electron density induced by the radiation field. The variation of electron density is 
a characteristic of the system, allowing one easily to connect classical and quantum 
descriptions of the excitation process, because charge density variation has the same 
meaning in both quantum and classical mechanics. A classical description of the electron 
density variation in a cluster is appropriate in the situation where plasmon excitations 
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dominate over the single-particle spectrum, because plasmon oscillations in clusters are 
an essentially classical effect. 

3.1. Basic equations 

Since our object of interest is the excitation of plasmons in metallic clusters, which 
have a distinctly classical nature, we now describe the collective motion of the electron 
density using Euler's equation and the equation of continuity. 

Euler's equation couples the acceleration of the electron density dv/dt with the total 
local electric field E acting on the density at the point (r, t) . It has the following form: 

^ = ^'> ™ 

The electric field E includes both the external field acting on the cluster and the 
polarization contribution arising from the variation of electron density. Expressing the 
total derivative on the left hand side of (0) as the sum of two contributions, arising 
from the change in velocity of the electron density in time and in space, one obtains: 



dt 



{v(r, t) ■ V}v(r, t) = --V<p(r, t) - -V / dr'^4f ( 17 ) 

m m J r— r 



Here <p(r,t) is the potential of the external field. The second term on the right 
hand side of ( |TTD describes the polarization force due to the variation of electron density 
8p{v,t). 

We assume that the external potential <p(r,t) is the solution of the wave equation. 
Therefore, we can put 

<p(r,t)=e iut <p(T), (18) 

where ip(r,t) satisfies the equation 

A<p(r) = -k 2 ip{r). (19) 

Here k = lu/c, c is the velocity of light, but in principle one can postulate a more 
complex dispersion law. We need consider only the positive frequency solution of the 
wave equation, because the formalism for the negative frequency solution is analogous 
to it. 

The total electron density in the cluster is equal to 

p(T,t)= Po (T) + 6p(T,t), (20) 

where p D ( r ) is the electron density distribution in a free cluster without an external 
field and 5p(r, t) is the variation of electron density caused by the external field and the 
polarization force acting together. 
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The motion of electron density in the cluster obeys the equation of continuity, which 
reads: 

dp(r,t) 



dt 



+ V-{p(r,t)v(r,t)} = 



(21) 



The simultaneous solution of equations (|TTD, ( p0|) and (^1|) with appropriate initial 
conditions and the initial distribution p (r) allow one to determine the variation of 
electron density 5p(r,t) as well as its velocity v(r,t). We solve this problem by using a 
perturbative approach on the external field ip(r, t). 



3.2. Perturbation theory 

It is easy to estimate the relative value of the first and the second terms on the left 
hand side of ([IT]). We see that the second term is negligible, provided the condition 
E -C mu 2 R/e is fulfilled. Substituting here the characteristic values to ~ 0.1, R ~ 10, 
one derives £<0.1 in atomic units or E <C 5 ■ 10 8 B /cm. Below, we assume that this 



condition is fulfilled and neglect the second term on the left hand side of (|17D, which 
means physically that the external field causes only a small spatial inhomogeneity in the 
electron density distribution within the cluster. In this limit, Euler's equation reduces 
to a Newtonian equation, which describes electronic motion in the cluster under the 
action of the external field and the polarization force. 

We express the solutions of ( |17D and (|2l|) in the following form: 



vnut 



(22) 



n=l 



v(r,f) = 5> n (r)e< 



(23) 



n=l 



By substituting these expansions into ([jj]), (0) and performing simple 
transformations, one derives 



v„ r 



— 5 ni iV^(r) + -^—X I r/r' 



mnuj 



mnu 



r — r 



n-l 



iujnSp n (r) + V • < p (r)v n (r) 



+ H V ' \ <W( r ) V n-fc'( r ) 
k'=l { 







(24) 



(25) 



Here 5 Ut i is the Kroneker symbol. One can exclude v n (r) from equation (^) by the 
substitution of (|24] ) in (|25|). Performing this transformation with the simultaneous use 
of ( |T9D and A|r — r'| _1 = — 47r5(r — r'), one derives the following equation: 

47re , A „ , , e_ ,Sp n (r') 



uon) 



\ e f 

■p (r) 5p n (r) + — Vp (r) • V / di 
J m J 



12 

e 

m 



5n,i (po{rMr)k 2 - V<p(r) • Vp (r)^ 

ft-i 

+iu ]T V- (^(r)v n _ fe /(r)) (26) 

k'=l 

The left hand side of equation (^) describes eigen-oscillations of the electron density. 
The electron density is almost constant within the cluster but varies rapidly near the 
cluster surface. Therefore, the terms proportional to p (r) and Vp (r) on the left hand 
side of (|26| ) determine the square of the frequency of the volume and surface plasmon 
oscillations respectively. The right hand side in ( ^6"D describes a driving force acting on 
the eigen- plasmon oscillations. 



The set of non-linear equations (|24|) and (|26|) must be solved iteratively. It is clear 
from the form of the equations that the index n corresponds to the order of perturbation 
theory on the external field <f(r). 

Indeed, for n = 1 equation (|2~6"D reduces to 
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III 



\ e f 

■ Po {t) 5p 1 (r) + -Vp (r)- V / dr' 



m 

and, for n = 2, one derives 



- (p (rMr)A; 2 - V^(r) • V Po (r)) (27) 



my m J 



= iwV-(«Jpi(r)v 1 (r)) (28) 

Equations (|27|) and (p8|) show that the variation pi(r') describes the linear response 
of the electron subsystem to the the external field <p(r), while P2(r') arises only in the 
second order of perturbation theory on <p(r), because pi(r') ~ <p(r) and vi(r) ~ <p( r )- 

Solving the set of equations (p4|) and (|26|) with <p(r) describing the dipole electron- 
photon interaction up to the n-th order, one can calculate the variation of electron 
density in the cluster caused by the field of n photons. 

The set of equations (^|) and ( ]2"6j ) is not confined in its application to photons. 
It can also be used to describe the dynamics of electron density under the action of 
any kind of external field, for example the electric field of a charged projectile colliding 
with the cluster. Indeed, by considering the partial spherical harmonic of the Fourier 
image of the Coulomb field of the projectile particle, one can derive from (27) the same 
expression for the variation of the electron density 5pi(r) as follows from the purely 
electro dynamical perturbative approach to the electron scattering problem [pit . 



3.3. Spherically symmetric case 

Equations (|2H|) and (p6|) are valid for an arbitrary shape of the initial distribution p (r). 
In the spherical case, the angular parts in ( |2~3| ) and ( |2"6| ) can be separated. Thus, the 
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cross section for n-photon absorption can in principle be extracted for arbitrarily large 
n, although the calculations become more and more tedious the higher n is. Let us 
consider this formalism in more detail. 

In the case of the spherically symmetric initial distribution, one can put p ( r ) = 
Po(r). This relationship allows one easily to exclude angular variables from equation 
(f26|). Using this relationship together with the partial expansion for 5p n (r) and (p(r), 

oo I 

6Pn(r) = E E WrWnr) ( 29 ) 

1=0 m=-l 
oo I 

¥( r ) = E E ¥i,m(r)Y^ m (n r ), (30) 

1=0 m =-l 



one derives 



((^) 2 " ^Po(r)) W, m (r) + ^0Y) I dr ' G i(ry)W, m (r') = 

e / \ e n_1 1 r 
= —S n ,i {po{r)^ m (r)k 2 - <p' ltrn {r)p' {r)) E ■ / ^n^^K) 



x jv<%(r) ■ V (^ Wil ^(r) + Jdr 



3=1 



-5 Pj {Y)(5 n - hl ^{r)k 2 + 4 7 rfy„_ i (r)) j (31) 

Here, Y/ i?Tl (n r ) is the sperical harmonic corresponding to the angular momentum I 
and the projection of the angular momentum m. When deriving (|31|), we have multiplied 
both sides of equation (p6|) by the sperical harmonic F z * m (n r ) and then integrated over 
dfl nr . We also used the well known expansion (see e.g. |27|j ) 



1 1 00 1 

-EE B^OWW) (32) 

-L 7 n ™ 7 



2Z 1 ■ , 



where function Bi(r,r') is defined as follows 

S,(r, r') = ^e(r' - r) + -^0(r - r') (33) 

Here 0(r' — r) is the step function. In equation (^l|) we have introduced the function 
Gi(r,r') = r 2 dBi(r,r')/dr, which is of the the form 

Gi(r, r') = l^&(r' - r) - (Z + l)^6(r - r') (34) 

When deriving (|3lf) , we have also made obvious transformations of the sum over k' on 
the right hand side of equation (p6|) using ([Tj^) and (0). Note that the sum over j 
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in ([Jl|) still contains the integrals over the angular variables. The integration over the 
angular variables in the sum is straightforward, but somewhat cumbersome. It is also 
clear that the non-gradient terms in the sum contain the integration of the product of 
three spherical harmonics, if one expands 6p n -j(r) and ip(r) according to ( |29| ) and (|30|), 
which we denote as 

ii(7,m|Zi,mi, \l 2 ,m 2 ) = J d^l nr Y l * m (n r )Y h , mi (n r )Yi 2tm2 (n r ) (35) 

The gradient terms in the sum contain the integration of a spherical harmonic 
multiplied by the scalar product of the two vector spherical harmonics: 



I 2 (l,m\li,m u \l 2 ,m 2 ) = \fh(h + l)l 2 (h + 1) / dVL nr x 

^(n^Y^^-Y^K). (36) 

This type of integral arises, when one expresses the gradients of the potential and 
the density according to (see e.g. P7|) 



V<H imi (r)Y"/ limi (n r 



= 5^ tmi (r)Y\-^M + y/h(h + lj-K^WY^K) (37) 

= <^(r) (n r ) + s]l 2 {l 2 + l)\<S>tJl (r) (n r ) (38) 

Here Yj m (n r ), Y^(n r ) are, respectively, the longitudinal and the transverse 
vector spherical harmonics, the definition of which one can find in |27|]. We mention 
some properties of these vector hamonics: Y^\n r ) = n r Y^ m (n r ), Y^(n r ) = 

VnnmW/fli + 1) and Yfcftn,) • Y« K) = 0. 



The potential _L (r) in (RSI) is as follows 



= + (2^1) J dr ' r ' 2B iM 5 Pt^ r ') ( 39 ) 

Since the explicit expressions for Ii(l, m\li, mi, \l 2 ,m 2 ) and I 2 (l, m\li, mi, \l 2 ,m 2 ) are 
somewhat lengthy, they are presented in |Appendix B| . 

Using the formulae written above, one can easily rewrite (^1|) in the following form 

((^) 2 " ^Po(r)) SptJr) + J dr'G^r^Jr') = 

-S n ,i (p (r)^, m (r)fc 2 - v'i, m (r)p' (r)) - 



e 
rn 



e ^ 1 



m 7^ n — 7 , 
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x { - 5p[ umi (r) (5 n - jtl <p h , m2 (r) k 2 + ^Sp^r))} - 
--^E— E i'2(i,m|Z 1 ,m 1 ,|i a ,m 2 )H, rai (0*aW (40) 

J=l «i> m l 

Here and below, we assume that the summation over h,mi and ^ m 2 is performed 
with the same limits as in (|29|) and (p0|). 

3.4- Surface and volume plasmons 

We now analyse equation (f40D and demonstrate that it describes both surface and volume 
plasmon oscillations. The surface and volume solutions of ( [40|) can be separated from 
each other, if one assumes that the initial distribution of electron density has the form 

Po(r) = ^Q(R v -r) (41) 

Here N is the total number of delocalized electrons in the cluster volume V = 4irRl/3. 

In this case it is natural to look for the solution of fl4Ci|) by expressing it in the 
following form 

*p£(r) = 5pff5(R s -r) + 5p v ^ (r)Q(R v - r) (42) 

In (|4lD and (f42|) we have introduced the two radii R v and R s and assumed that 
R v < R s = R, but R v — i> R s = R, where R is the cluster radius. Such a relationship is 
necessary for the elimination of the uncertainties, which arise in ( f42|) and the subsequent 
formulae in the vicinity of the cluster radius. 
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Substituting ([41]) and fl4"2|) in ([4(]), performing straighforward but lengthy 
calculations of the integrals, one derives the following equation 

((urn) 2 - u>}) 5p£ ) (r)e( J R-r) + 

+ ((uju) 2 - ujf) 5p s ^5(R - r) = 

_ 47reiV(/ + l) f R dr'r n+2 6o v{n) (r') + 

Ne 2 Ne 2 
+ 5 n . mm (r)k 2 Q(R -r) + —6 n ,M m (R)S(R - r) + 
mv mv 



e „,„ ^ 1 



x 



J = l il,m! 



+5 n _ J -iA; 2 5pf r i i ^ 2im2 ( J R)| - 
e " _1 1 

e ™ _1 1 

~mJ¥^ R ~ T ^^ ~^ZT- ^ /2 ^' m ' Zl ' mi ' m2 ) X 

^ j = l ^ 3 li,mi 

x^fil (^^(i?) + j^-j f* dr'r*B l2 (Ry)6 P ^\r')\ - 

e n_1 1 
9 - r) V ; V I 2 (l, m\h, mi, |/ 2 , m 2 ) x 

J = l J l 1 ,m 1 

h,m,2 



When deriving we have used the fact that Q'(R — r) = —5(R — r). Also we 
have introduced the function 



^'2—2 ^/l2+2 



Di(r, r') = l 2 (l 2 - l)- s ^ T e(r' - r) + (Z 2 + 1)(Z 2 + 2)-^9(r - r') (44) 



17 



The left hand side of equation (^) describes volume and surface eigen-oscillations 
of the electron density characterised by the angular momentum I. The surface plasmon 
resonance frequency uj\ is the same as in (0). The volume plasmon resonance frequency 
is equal to 



^ = i^r (45) 



In equation fl43|), u p appears in expressions involving I, but is nevertheless is independent 
of I, as one sees in equation (f45|) . The physical reason for this is that the volume plasmon 
oscillation is degenerate with I. The right hand side of fl43|) describes a driving force 
acting on the eigen-plasmon oscillations. 

Surface and volume terms on the right hand side of equation (|43|) have not been 
regrouped, in order to stress their correspondencce with terms in (|40|). It is seen from 
([43|) that equations for the volume and surface plasmon oscillations can be separated 
and will then read as follows: 

e n_1 1 

V : V I 1 {l,m\l 1 ,m 1 , \l 2 ,m 2 ) x 

j = l J l 1 ,m 1 

e n_1 1 

p Y\ : Y\ h(l, m\h, mi, \l 2 , m 2 ) x 

mr z ~i n — ] , 



4ireN(l + l) f R , /l+2 v(n) 



m(2l + l)VR l + 2 Jo 
e ^ 1 



H Y] : V] h(l,m\h,mi, \l 2 ,m 2 ) x 

m 7^ n - i , 



18 



+5 n ^ j:1 k 2 5p s l { i 3 X i (p him2 (R)^ - 



e ^ 1 



~^R2 ~^ZT- ^ /2 ^' m ^ mi ' l^ 2 ' m ^ x 

^n-i.i^W + T^T) /"d^B^r^S^w) (47) 



The set of non-linear equations ( |4q) and fl47|) must be solved iteratively starting 
from n = 1. It is clear from the form of these equations that the index n corresponds to 
the order of perturbation theory on the external field <p{r). It is also seen that the chain 
of equations for the volume plasmon oscillations is self-sufficient, while the solutions 
for the surface oscillations also depend on those for volume, which is physically clear, 
because the electric field caused by volume oscillations of electron density must influence 
surface oscillations. 

From (|46|) and (0), one derives the following solutions for n = 1: 

S f) {r) = ^fcVUr) (48) 



2 



These equations show that, if the external field (|30| ) is characterized by a certain 
angular momentum l a and its projection m , i.e. tpi, m {r) = 5i y i 5 m , mo ipi , mo (r), 
then S PiS = 6 Uo$m,mJp s £L and Spfflir) = S ljlo S mimo Sp v ^l o (r). Assuming these 
dependencies, the solutions of ( [46]) and (|47|) for n = 2 read as: 



$Pi,m ( r ) = TTTT"^ osh(l,m\l ,m , \l ,m ) x 

x |^S' (r) ^ o , mo (r) + / o ^r'GJr,r')<l(r0) - 



*<W£» (vi a ,m (r) + £ ! ^B io (r,r')^(r')j (50) 
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jj'P) - AireN V ± ^ f R dr'r' l+2 6o v{1) (r') + 



o 



-Ji(/,m|/ ,m , \l ,m ) x 



m((2cu) 2 - cj 2 ) 

+fc 2 5p / s « o ^ o , mo ( J R)}- 

-I 2 (l,m%,m , \l ,m ) x 



m J R 2 ((2cu) 2 -w/ 



x*pfiL ^,m.(i2) + [ dr'r*B l XRy)5 P i%y)j (51) 

By performing similar transformations, one can find the solutions SpJ^ (r) and 
for arbitrarily large n, although the formulae become more and more tedious the larger 
n becomes. These formulae demonstrate that, in the higher orders of perturbation 
theory, plasmon resonances with angular momenta larger than the angular momentum 
of the external field can be excited. Indeed, the selection rules for the integrals 
m\l , m , \l ,m ) and J 2 (/, m|/ , m Q , \l ,m ) (see |Appendix B|) show that the angular 



momentum in 8p^ (r) and can be twice as large as l a . Equations (|46| - [51]) also 

demonstrate that the plasmon resonances in 8p^\r) and 5p s l ^ arise when to = uj p /n 
and uj = tui/n respectively. These equations indicate a significant shift of the plasmon 
resonance profiles towards lower frequencies in the highest orders of perturbation theory. 

These results have a simple physical explanation. Absorption of several quanta of 
the external field (photons) by the cluster leads to the excitation of non-dipole plasmon 
oscillations of the electron density. 



3.5. Fast electron- cluster collisions 

Equations ([46]- |5l]) can be used for the analysis of the balance between the surface and 
volume plasmon oscillations in the cluster. We demonstrate this for the example of fast 
electron scattering on a metal cluster. In this case, the external field of the projectile 
electron can be characterized by the Fourier component of the Coulomb potential 

<p{r) = %e^, (52) 

where q = p — p' is the transferred momentum of the scattered electron. 
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The partial expansion of this potential reads as: 

oo m=l 

^(r) = 4tt £ E i l fi,m(r)Y* m (n q ), (53) 

i=0 m=- 

where the partial component of the potential (pi >m (r) is equal to: 

47T 

<Pi,m(r) = ^-ji(gr)F Z!m (n r ), (54) 



and ji(qr) is the spherical Bessel function (for definition see e.g. 

The form and properties of (pi >m (r) - (0) are exactly the same as assumed in (O) 
and Q3"0|). Therefore, from and (ff9f), one can immediately derive: 



«£?<r) = ^ (55) 



and 

= ?R^F) i J ' ;(gi?) " (^-^)(2i + i) X dr V J/(gr ' ) | (56) 

In the case of inelastic electron scattering, u has the meaning of the transferred 
energy in the collision, u = Ae — e — e'. Calculating the integral in fl55| ) with the use 
of the well known properties of spherical Bessel functions (see e.g. [f25j), one derives 

s(1) {2l + l)u:f 3l {qR) u% . 
^ = w2 _ wj! g2i? - 1 j^^m (57) 

Expresssions ( j55f) and (|57|) coincide with those calculated in [^j], in the plasmon 



resonance approximation by purely electrodynamic means, as the response of a dielectric 
sphere, having dielectric permeability e = 1 — 9(R — r)uj 2 /u^. 

From (B^) and (0), one can easily elaborate the electron inelastic scattering cross 
section in the plasmon resonance approximation, using the method described in [21|: 

V* £(2i + lft , faii) ^!L_ + (58) 



cfe'dfl vrpg 4 V ~ ^ 2 ) 2 + Ae 2 rf 2 



vrpg 2 Y - ^ P 2 ) 2 + Ae 2 Tf 



x 



x - ji+i{qR)ji-!{qR) - ^j l+1 {qR)ji{qR)^ 

This cross section is totally determined by collective electron excitations in the 
cluster. The first and the second terms in (|58|) describe the contributions of the surface 
and the volume plasmon excitations respectively. In ([58]), we have also indroduced 
widths, rf and r", of the surface and volume plasmon resonances. They originate from 
the Landau damping mechanism of the plasmon excitations. For their determination 
we refer to the recent paper [BTf] . 
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3.6. Multiphoton absorption 

Next, we apply equations P6| ) and (fj7|) to the description of the multiphoton absorption 
process. In this paper, we focus our consideration on the analysis of plasmon excitations. 
If surface or volume plasmon resonances are excited by photons, i.e. u ~ u p , then it is 
easy to check that the following condition is fulfilled uR/c ~ u p R/c <C 1, where c is the 
velocity of light. This condition implies the validity of the dipole approximation. 

In the dipole approximation, one can neglect the momentum of the photon and put 
k = 0. In this case, equations ([R|) and (|47|) are simplified dramaticaly. Indeed, from 
(fSf), one derives 

8p?2S (r) = (59) 

s(n) 



This result also simplifies equation (37) significantly for Sp^ . After some trivial 
transformations it reduces to: 

p n_1 ft 



„ i n — 7 , 



m 

.; ■ ' 

e 2=J 5„ 



_j_ *LJ>L ^2 J 2 (Z,m|Zi,mi, l^^fyS^^OR) ( 60 ) 
miir r— ; n — 7 , 



h,m,2 



The partial component of the linearely polarized dipole photon field is equal to 

/47T 

fi,m{r) = -]l —Er5i t i5 mfi (61) 



Here E = JlitTno /V is the strength of the photon's electric field and V a is the 
normalization volume of the photon mode. Substituting (|5"T| ) into (fT7|), one derives 

/V ^ 2\ r s(n) [^ Ne 2 E 

{(ujn) - w, J bp^J = -W y^^5„,id/,id m ,o + 

_ (62) 

47rejB £ / 2 (i,m|i ll m ll |l,0)H ( ^ 1) (63) 

This equation should be solved iteratively starting from n — 1. For n = 1, the single 



non-trivial solution, Spi , reads as 
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Then, for n — 2, the solution of ( p^ ) is of the form 

^ = i/f mfi (S' < - M ?) wm|i - o ' i - o > (65 » 

The selection rules for ^(l, m\l, 0|1, 0) (see |Appendix B|) show that this integral 



does not vanish, when / = and m = or I = 2 and m = 0. Therefore, for n = 3 from 
(|B3D, one derives 

r s(3) = /4tt eE 
Pl ' m V 3 mR{{3uj) 2 - uj?) 

{h(l, m|0, 0|1, 0)5pJ 2) + J 2 (/, m|2, 0|1, 0)5pg } } (66) 

In ( jB6|) , only the second term in brackets gives a non-zero contribution, since 
J 2 (i, m|0, 1 1 , 0) = (see |Appendix B| ). Substituting 5poo an d <^°20 fr° m (|S5|) and 



^Pi'o'' from (|64D in (|66|) and using the explicit expressions for the angular integrals given 
in [Appendix B , one obtains 

K l/2 Ne 3 E 2 



s(2) 



<Wo 



X s ( 2 ) 



X s ( 3 ) 



3m 2 RV uj 2 (uj 2 -ujf) 
An 1 / 2 Ne 3 E 2 



3V5m 2 RV (uj 2 - u 2 )((2uj) 2 - u 2 ) 
64tt\ 1 / 2 16 + 3^ Ne 4 E 3 



3 ) 75m 3 R 2 V {uj 2 - uj 2 ){{2u) 2 - cu 2 2 )((3cu) 2 - uj 2 ) 



Sp »&> = 1/2 12(2 + 75) Ne*E* 

P3 '° \ 7 J 75m 3 W {uj 2 - oj 2 )((2lu) 2 - uj 2 2 )((3uj) 2 - uj 2 ) 1 ' 

4. Induced multipole moments in the cluster 

Let us now calculate the multipole moments of the cluster induced by an external 
radiation field on the basis of the model developed in the previous section, and analyse 
their plasmon resonance structure. 

The induced multipole moment of the cluster is equal to 



Qi, m = y / dVr% m (n r )5p(r) (68) 

where the variation of electron density Sp(r) is determined in (|2~9|) and (fj2]). Subsituting 
(f2~9|) and ( |42"D in (|68|) and putting Sp"^ (r) = for any n in the dipole approximation 
as follows from (|59[), one derives 
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Substituting in d69|) Sp^Q from fl64|), one obtains the expression for the dipole 
moment of the cluster induced in the single-photon absorption process 

Df)( M ) = qW = - 2 7Ve ^. (70) 

The explicit expressions for the partial electron density variations 5pf^ entering 



(|69|) for n = 2 (two photon case) and n = 3 (three photon case) are given in (|67|) . 
Subsituting the partial electron density variations Spo^ and ^p^cf from ( |67|) into (|69|) , 
one derives the expression for the monopole and quadrupole moment of the cluster 
induced in the two photon regime 



''° 2m 2 R 2 uo 2 {uo 2 -ul + \DT t ) 

nV){ \ — - —- Ne 3 E 2 

W [U)-V2,o- 5m 2( u 2_uf + iuJ r 1 )((2u) 2 -u 2 2 +i2u;r 2 ) 1 } 

Here, we have introduced the plasmon resonance widths Ti and T 2 which take into 
account Landau damping of the dipole and quadrupole surface plasmon resonances. 
They must be determined separately, e.g. by an ab initio approach (see |2T| ). 

By absorbing three photons one can induce dipole and octupole moments in the 
cluster. Substituting Sp{^ from p7| ) into (|55|), one derives the expression for the induced 
dipole moment 

X {uo 2 - uo 2 + iuoT 1 ){{2uo) 2 - J 2 + i2uoT 2 ){{3uo) 2 - uo\ + i3^rj 72) 
The expression for the octupole moment induced by 3 photons following from (|37j) 
and (|69|) reads as 



n(3) L , n(3)_ 12(2 + ^5) 



X {uo 2 -uo\ + luoT^uo) 2 - uo 2 + i2uoT 2 ){{3uo) 2 - uo 2 + i3^r 3 f ?3) 
Here, we have also introduced the octupole plasmon resonance width T 3 . 
Expressions (|70|-[f3|) demonstrate that the multipole moments induced in the 

cluster during multiphoton absorption processes possess a prominent plasmon resonance 

structure. The nature of these resonances is the same as occurs in the mutiphoton 

absorption cross sections discussed in section |2|. 

The connection between from (|70|) and the cross section a\ found in (|J) is 

straightforward 

ai = ^ ImD W^) (74) 
cE 
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In the multiphoton regime, the connection between the induced multipole moments 
of the cluster and the multiphoton absorption cross section becomes more complex, 
which is apparent from the the classical nature of the expressions (f71]-[73|) and the explicit 
dependence of the multiphoton absorption cross sections on Planck's constant. The 
discussion of this interesting relationship is however beyond the scope of the present 
paper. 

5. Conclusion 

In this paper, we have developed a formalism which allows one to calculate the cross 
section for multiphoton absorption in the plasmon resonance approximation. We 
have demonstrated that plasmon excitations with angular momenta larger than 1 
substantially alter the profile for multiphoton absorption as compared to the single- 
photon case. 

Our model is formulated in terms of a charge density distribution function p(r) 
for the cluster. This means that, in principle, one can study the response for different 
charge density profiles including deformed ones. Our model is a semi-classical one, which 
neglects the granularity of charge in the system. This is consistent with the principles 
underlying the jellium picture. It is appropriate for metallic clusters and, to a lesser 
extent, for fullerenes. 

In the classical formulation of our model, we have used Euler's equation for 
hydrodynamic flow, together with the equation of continuity. We have demonstrated 
that the results following from our model are consistent with direct estimates of the 
matrix elements for the multiphoton absorption process. The theoretical formalism we 
have developed is not confined in its application to photons. It can also be used to 
describe any kind of higher order plasmon excitation processes, for example those which 
arise by multiple scattering of electrons within a cluster. 
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Appendix A. Matrix elements of plasmon resonance transitions 

In this appendix, we evaluate the matrix elements of plasmon resonance transitions in 
the plasmon resonance approximation by the use of the sum rule. 

For a stationary (i.e. explicitly independent of time) operator F of an observable 
physical quantity, characterising a system of particles with the Hamiltonian H, one can 



25 



formulate the following sum rule (see e.g. p5*|) 



]T^o|(n|F|0>| 2 = i(0| 



10} 



(Al) 



Here u n o = e n — eo, the summation is performed over all excited states of the system 



and H, F denotes the commutator of the operators H and F: 
Applying the sum rule ( |AT| ) to operator F defined as 

F = £F(r fc ), 



H, F 



HF - FH. 



(A2) 



where F(rk) is a function of the coordinates of the k — th electron and the summation 
in (|A2|) is performed over all particles in the system, one derives 



$> n0 |(n|F|0>| 2 = ^ / rfr|VF(r)| 2 p(r) 



(A3) 



Here p(r) is the ground state charge density distribution in the system. Applying 
the general rule (|A3|) to the function 



F(r) 



4tt 



y 21 + 1 

and performing the integration in 



r l Y lm {n) (A4) 
with the density distribution (^1|), one derives 



J^^nO (n\Ql m \0) 



-u?R 2l+1 



(A5) 



Here, we have used the following definition of the multipole moments operator 



47T 



N 



Qim = y Yl e*r^ m (n k ). 



(A6) 



The plasmon resonance frequencies in ( |A5| ) are defined according to (0). 

Using (|A5|), one can easily evaluate the matrix elements of plasmon resonance 
transitions in the plasmon resonance approximation. Indeed, assuming that plasmon 
excitations dominate in the sum over n in (|A5|) , one derives 



Qim = (n\Qim\0) 



lhu t R 21 - 



eR 



i-i 



Nh 
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2mui V 21 + 1 



(A7) 



Equation ( |A7| ) gives the matrix elements of plasmon resonance transitions for 
an arbitraty large angular momentum I. The correctness of the result (|A7|) can be 
independently varified by performing calculations of the multipole dynamic polarizabilty 
of the cluster in the plasmon resonance approximation. Indeed, using (|A7|), one derives 
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which is the known expression for the dynamic multipole plasmon polarizability in the 
plasmon resonance approximation (see e.g. |T0, T7|). 

For the dipole plasmon resonance transition, one derives from ( |A7| ) 



^ Nh ,. , 

Qw = d w = ez w = eW (A9) 

V Imuji 

which is consistent with the dipole sum rule (^). 

Equation ( |A7|) allows one to evaluate the matrix elements of electronic transitions 
between various plasmon resonance states. To demonstrate this, let us rewrite ( |^7| ) in 
the form 



Qlm = \ 7^-T I drr l+2 p l0 (r) (A10) 



2/ + 1 

Here, we have introduced the radial transition density pzo( r ) between the ground 
state and the excited plasmon resonance state with angular momentum I and used the 
relationship I\(l n , m n \l, m|0, 0) = o~i ny ib~ mnym l 'y/An, when calculating angular intergrals in 

The radial transition density pio(r) is localised in the vicinity of the cluster 
surface. Qualitatively, this is clear because pm(r) describes the plasmon excitation. 
Quantitatively, this was proved by ab initio computations of the transition densities in 
the Na^o and Nag2 clusters within the jellium model in [0. Therefore, to a reasonable 
accuracy, one can approximate pio(r) by the delta function 

Pio(r) = Pio5(r - R) (All) 



Substituting (fATTD into ( \KWj ) and comparing the result of the calculation with (|A7|) , 
one can determine the value pio entering ([All ). The result of this calculation reads as 



e 31HN /A s 

"» = 7pfc (A12) 

Let us now evaluate the matrix element for the dipole transition between plasmon 
resonance modes, which reads as 




z hh = \hr drY i,o(n)rpi 2 h(r) (A13) 



Here, pi 2 ^ (r) is the electron transition density between the dipole and the quadrupole 
plasmon modes. This transition density can be evaluated via the transition densities 
Ph,h{ r )i Pho( r ) an d the ground state electron density of the cluster p 00 = eN/V as 
follows 

Poo 
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Substituting ( |A14| ) into ( |A13| ) and performing simple transformations, one derives 



Uirhh f h(2l 2 + l) \ 1/ \ ^ h xfn s (a,k-\ 

Zhh = e\— — — — — — 4r\ h{l 2 ,m 2 \l^\l ll m 1 ) 8(0) A15 

V 3 \h{2h + 1)/ 

This equation has an uncertainty, which originates from the fact that we have 
assumed zero width for the domain in the vicinity of the cluster surface in which plasmon 
excitations take place. By introducing a finite width AR for this domain and using 



one of the standard representations of the 5-function |25j to resolve the uncertainty, 
8(0) ~ 2/irAR, one finaly derives 

^ = 8e V^(M277Ti)J /l(/2 ' m2|1 ' 0|/l ' mi) ^A^ (A16) 

The explicit expression for the angular intergral h(l 2 , m 2 \l, 0|/i, mi) is given in 
Appendix B| . In the case of the transition between the dipole and quadrupole plasmon 



resonance states this integral is equal to I±(2, 0|1, 0|1, 0) = — l/v57r. Substituting this 
value into ( |A16|) and performing simple algebraic transformations, one arrives at the 
expression for the matrix element describing the transition between the dipole and the 
quadrupole plasmon resonance modes. 

^ = - e 3U) ^AR (A17) 
Appendix B. Integrals Ii(l,m\li,mi,\l 2 ,m 2 ) and I 2 (l,m\li,mi,\l 2 ,m 2 ) 

The angular integral, 

I 1 (l,m\l 1 ,m 1 \l 2 ,m 2 ) = J dtt nr Y^ m (n r )Y lltmi (n r .)Y l2 , m2 (n r ), (B18) 
is well known and can be found in many textbooks (see e.g. [pTfl). It is equal to 



^i 1+ i 2 -i.im + mh + i)(2h + i: 



h^mlh.m^m,) = (-l)"V^- y il i_ l± x 

x / I h l 2 \f I h l 2 \ 
\-m Trt\ m 2 J \ J 

Here, the integral is expressed via 3j-symbols (for definition see e.g. ||27|| ). 
The angular integral, 

l2(l,m\h,mi\l 2 ,m 2 ) = ^h(h + l)l 2 (l 2 + 1) / dfi nr ,^* m (n r ) x 



(B19) 
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can be expressed via the sum of products of 3j-symbols and 6j-symbols (for definition 
see e.g. |^7||) after performing the following transformations. 

Using the standard relationships for spherical vector harmonics, written in [27| on 
page 210, one derives 

I 2 (l,m\l umi \l 2l m 2 ) = yjh(h + l)l 2 (l2 + l) J dn n Y* m (n r ) x 



x 



/ h + 1 - 
2k + 1 

1 h + 1 - 



Wi-l 



n r ) + 



_J^_Y l2+1 
2l 2 + l -h,m 2 K"rJ-r y 2/2 + 1 h,mx 



n r ) + 



2 h + i'h, mi ^r 



(B21) 



The integrations arising in (|B21| ) can be performed and expressed via the sum of 
products of 3j-symbols and 6j-symbols, using the standard formulae given in f27j on 
pages 222 and 236. The result of the calculations of these integrals reads as 



l2(l,m\li,mi\l 2 , m 2 ) = 

= (-l) h+h+l+m+1 ^h(h + 1)Z 2 (Z 2 + 1)(2Z + 1) 



mi m 2 —m 



x 



x • 



, (/ 1 + l)(/ 2 + l)(2/ 1 -l)(2/ 2 -l) I 2 -l Afli-ll 2 -ll 

4vr V i 1 l 2 h 1 



/ / 1 (/ 2 + l)(2/ 1 + 3)(2/ 2 -l) (h + 1 l 2 -l I fuiij-nl 
V 4tt J I l 2 h 1 j + 



/ (/ 1 + l)/ 2 (2/ 1 -l)(2/ 2 + 3) / 2 + l I a-1 Z 2 + l n 

V 4tt i I Zi Z 2 1 J 



+ 1 



/ /i/ 2 (2Zi + 3)(2Z 2 + 3) / h + 1 Z 2 + 1 Z \ f Zx + 1 Z 2 + 1 Z 
4vr V ) \ l 2 hi 

For the particular cases of interest, one derives from ( |B22| ) 
J 2 (1,0|0,0|1,0) = 
7 2 (0,0|1,0 



(B22) 



/ 2 (2,0 

/ 2 (i,o 

I 2 (3,0 



1,0 
2,0 
2,0 



1,0) 
1,0) 
1,0) 
1,0) 



1 












1 












iov? i 1 + 


16 




3^ 


) 


18 


(■♦ 




5V105tt 


2 
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